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ABSTRACT 

We present a state-of-the-art scenario for newly born magnetars as strong sources 
of Gravitational Waves (GWs)in the early days after formation. We address several 
aspects of the astrophysics of rapidly rotating, ultramagnetized neutron stars (NSs), 
including early cooling before transition to superfluidity, the effects of the magnetic 
field on the equilibrium shape of NSs, the internal dynamical state of a fully degenerate, 
oblique rotator and the strength of the electromagnetic torque on the newly born NS. 
We show that our scenario is consistent with recent studies of SNRs surrounding 
AXPs and SGRs in the Galaxy t hat constrain the elect romagnetic energy input from 
the central NS to be 10^^ erg ()Vink fc Kuipeil 120061) . We further show that if this 
condition is met, then the GW signal from such sources is potentially detectable with 
the forthcoming generation of GW detectors up to Virgo cluster distances where an 
event rate ~ l/yr can be estimated (Stella et al. 2005). Finally, we point out that the 
decay of an internal magnetic field in the 10^^ G range couples strongly to the NS 
cooling at very early stages, thus significantly slowing down both processes: the field 
can remain this strong for at least 10'^ yrs, during which the core temperature stays 
higher than several xlO^ K. 
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1 INTRODUCTION 

Gravitational wave emission from spheroidal, rapidly rotat- 
, ing, isolated N Ss has long been conside red in the astrophys- 
ical literature (jOstriker fc Gunnlll969l ). A natural origin of 
this distortion can be the anisotropic pressure from the in- 
ternal magnetic field. However, before the early nineties, the 
inferred magnetic fields of NSs were ~ 10^^ G, which implied 
tiny deviations from spherical symmetry and, thus, weak 
GW emission that would be detectable, in principle, only 
from the nearest sources and wi th years-long observations 
l|Bonazzola fc Gourgoulhonll 19961 ). 

The discovery of the Soft Gamma Ray Repe aters (SGRs) 
and Anomalous X-ray Pulsars (AXPs) (cfr. iMazets et al.l 

I1979I . iMereghetti fc Stellalll995l ) led to the idea that these 
peculiar high-energy sources could be ultramagnetized NSs, 
magnetars, with external (dipole) fields in the 10^* -r- 10^^ G 
range an d with internal fie lds at least one order of magnitude 
stronger (|Paczvnskvll992l . Duncan fc Thompson 1992, DT92 
hereaf t er, Th o mpson fc Duncan 1993, 1995, 1996, hereafter 

lTD93l : lTD95l : ITD96| . respectively). The core of the proto- 
neutron star (PNS) experiences a phase of neutrino-driven 
turbulent convection (during a few tens of seconds). Strong 
differential rotation is also present if the nascent NS is spin- 
ning at millisecond period (we indicate the initial spin, in 



milliseconds, as Pi/ms), providing a total free energy reser- 
voir of up to ~ 10^'^(Pi/ms)~^ erg. The combination of 
these two factors can power an a — Q dynamo that acts 
coherently over the whole NS core l|DT92l ). A seed magnetic 
field is twisted into a mainly toroidal configuration and its 
intensity amplified to values ~ lO^'' G. The total energy 
available from the core differen tial rotation correspo nds to 
a maximum field of ~ 10" G (|TD93l : iDuncanI [l993 ) . Mil- 
lisecond spin periods are required for differential rotation to 
play an important role and to achieve magnetic field coher- 
ence lengthscales comparable to the stellar radius. In more 
slowly rotating NSs differential rotation is almost negligible 
and a convection-driven, a-type dynamo results that acts 
stochastically and leads to a much weaker large scale field 
(|DT92l : lTD9a ). In this formation scenario, magnetars repre- 
sent those NSs that are born very rapidly rotating, at 5C 3 
ms spin period. 

As the magnetar model received increasing support from 
observations during the last decade it was shortly after real- 
ized that such objects at birth represent promising sources of 
GWs, be cause of their lar go magnetic deformation and rapid 
rotation (|lokall200ll : [Paio mba 2001). These early suggestions 
assumed that a simple dipole field extended throughout the 
NS and its magnetosphere thus causing NS deformation and 
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magnetodipole spindown. ICutleil (|2002l ) first highlighted the 
crucial role of the internal field structure in the GW emis- 
sion efficiency of a magnetically distorted NS. He pointed 
out that a NS with a strong internal toroidal field has a pro- 
late shape and would thus provide the best chance of being 
a strong GW emitter. Conversely, GW emission from a NS 
with an oblate distortion (such as that caused by a dipole 
field) would be rapidly quenched, since such an object would 
tend to spin around its symmetry axis (see §[21)- Based on the 
energetics and likely recurrence time of Giant Fla r es suc h as 
the Dec 27 event from SGR 1806-20, IStella et alj (|2005l ) de- 
rived a lower limit of ~ 8 x 10^^ G for the strength of the 
core magnetic field in this object. This suggests that internal 
fields of magnetars may be even str onger than previousl y 
thought, rea c hing the lO '^^ G range. [ Lazzati et al.l l|2005l ). 
iNakar et al.l (|2006l ') and IPopov fc SternI l|2006t ). based on 
the apparent lack of Giant Flare-like events in the BATSE 
archive, suggest longer recurrence times for such events, 
thus reducing the m inimum energy requirement proposed by 
IStella et al.l (|2005l ). However, the correlation betwee n short 
GRBs and local galaxies found bv lTanvir et al. I ('2005') leaves 
room for a fraction of such events to be associated with 
a pop ulation of young NSs. Furthermore, iKaminker et al.l 
l|2007l ) showed that, if the thermal emission observed from 
AXPs is indeed powered by the decay of a magnetar's B- 
field, then the field strength in the NS crust must be ~ 10^® 
G. 

Fields ~ 10^^ G would imply remarkably large (prolate) 
magnetic deformations, es ~ 10~^ (see § I2.1.2|l which, in 
turn, would make GW emission at birth from such objects 
strong enough to be detectable up to the Virgo cluster dis- 
tance, with Advanced LIGO/ Virgo-class detectors. The in- 
tegrated mag netar formation rat e in Virgo is estimated to be 
1 per year (IStella et al.l 



inforcing the case for fast spins at birth. 
If magnetars are born with millisecond spin periods (here- 
after millisecond magnetar), their initial spin energy will 
be £;,pi„ ~ 3 X 10^^(Pi/ms)-^ erg. Here and throughout 
we assume a typical radius of 12 km and a mass of 1.4 
M0. The spindown timescale through magnetic dipole ra- 
diation will be very short, ~ one day for external dipole 
fields in the lO" G range. Most of the initial spin energy 
would thus be rapidly transferred to the surrounding su- 
pernova ejecta (| Allen fc Horvathll2004l '). Recent studies have 
also shown that, for newly formed magnetars with dipole 
fields in excess of ~ (6 -^ 7) x 10^^ G and spin period of « 
1 ms, strongly magnetized, relativistic winds are produced 
which can carry away most of the initial spin energy in a 
matter of minutes, thus opening the possibility of an even 
faster transfer of energy to the ejecta ("Thom pson et all2004l : 
[Bucciantini ct al. 2006; Motzgcr ct al. 2007). In both cases, 
present-day SNRs around known magnetars should bear the 
signature of this larger-than-usual energy injection. For ini- 
tial spin periods less than 3 ms, the injected energy would 
be > 3 x 10^^ erg, making these remnants significantly more 
energetic than those surrounding ordinary NSs. 
The X-ray spectra of the SNRs surrounding known magnetar 
ca ndidates (two APXs an d two SGRs) have been analyzed 
bv lVinkfcKui^ iooi) who found that the total energy 
content of the ejecta in these remnants does not appear to 
be different from the energy in remnants surrounding com- 
mon NSs (« 10^^ erg). This result implies that either these 
NSs were not born rapidly rotating, thus challenging the 
a — n dynamo scenario, or their initial spin energy must 
have been lost without a ppreciably energizing the surround- 
ing ejecta. iDall'Osso fc Stella (2007i ) discussed the possibil- 
ity that most of the initial spin energy of millisecond magne- 
^ ^ . ^^^^ is realeased through gravitational radiation. They con- 
year in the Galaxy alone ( Gaensler et al.|[l999t ). thus mak- el uded that, in order to account for the constraints derived 

bv lVink fc KuipeJ (|2006l ) , the internal toroidal field at birth 
should be in the 10^*^ G range, consistent with the inference 
by IStella et al.1 Hooi) for the field strength in SGR 1806- 
20 and, by extension, magnetars in general. A similar con- 
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ing newly formed magnetars potentially interesting sources 
for next generation GW detectors, with a fairly high rate of 
occurence. 

As an alternative scenario, it was proposed that strongly 
magnetic NSs can form as a direct result of magnetic flux 
conservation during core-collap se of a massive star with un- 
usual ly s trong magn etic field (jPerrario fc Wickramasinghd 
l2006l . cf. IUsovlll992l ). In this picture, magnetars need not 
be born rapidly rotating: their magnetic field is stronger 
than in ordinary NSs because the magnetic dipole mo- 
ment of their progenitor stars was particularly large. 
iFerrario fc Wickramasinghd (|2006l ) have shown that a pop- 
ulation of NSs with dipole fields up to ~ 10^^ G can be 
obtained in this framework starting from a realistic Galac- 
tic population of high-mass, high-field stars. We note, how- 
ever, that the upper end of their field distribution seems 
to fall short of the minimal requirement within the mag- 
netar model. Internal fields of at least several times 10^^ 
G are needed bot h to ex plain the overall energy output of 
SGR s and AXPs (|TD95| . Thompson fc Duncan 2001, here- 
after ItDoJ), and to make the magnetic field decay timescale 
comparable to, or sh orter th an, the estimated ages (~ 10^ 
yrs) of these sources |TD96|). 

iGeppert fc Rheinhardtl (|2006l ) have shown that a field am- 
plified to magnetar-strength in a newly formed NS will sur- 
vive an early unstable phase (and subsequent dissipation) 
only if the NS spin period is shorter than ~ 5 ms, thus re- 
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iBucciantini et al] (|2007l ) have also studied the possibility 
that the huge spin energy of a millisecond magnetar is 
promptly released in the form of a highly coUimated wind, if 
the external (dipole) field exceeds 10^^ G. The wind breaks 
through the surrounding supernova ejecta and produces a 
GRB-like event without directly energizing the ejecta. How- 
ever, the currently estimated magnetar formation rate ex- 
ceeds the rate of observed GRBs by 2 orders of magnitude. 
Hence, this process can i nvolve at most a tiny fraction of 
newly formed magnetars (|Bucciantini et al]|2007l ). 
Thus, a number of astrophysical arguments point to possi- 
ble field strengths in the 10^^ G range in magnetar interiors 
and suggest a possibly relevant role of GW emission in the 
early evolution of these objects. The aim of this paper is 
to explore the physical consequences of this hypothesis and 
draw a state-of-the-art scenario, in light of all approxima- 
tions and uncertainties. The paper is organised as follows: 
in § [2] we introduce the model for the early spindown of a 
magnetar subject to both magnetodipole and GW torques 
and study the main factors determining the GW emission 
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efficiency of a newly formed magnetar. In § [3] we calculate 
the expected energy loss through GW emission from newly 
born magnetars as a function of the relevant NS parameters 
and re-ad dress the problem of signal detection already dis- 
cussed by IStella et al.l l|2005h . Finally, in § [5] we address the 
problem of the thermal and magnetic evolution of NSs, ex- 
tending previous treatments of field decay to field strengths 
~ 10^^ G. 



2 MAGNETAR EARLY SPINDOWN: 

ELECTROMAGNETIC VS. GW EMISSION 

A newly formed NS with a strong toroidal magnetic field 
has a prolate deformation induced by the field (es, de- 
fined in § I2.1.2|l with symmetry axis along the magnetic 
axis. If the magnetic and spin axes are misaligned, with tilt 
angle x, the angular velocity vector (ft) precesses (in the 
NS frame) arou n d the fixed angu lar momentum vector (L). 
ICutler fc Jonej (|200ll ) (see also ICutled |2002| ) studied the 
precession dynamics taking also into account the effects of 
the centrifugal deformation of the NS. They showed that the 
triaxial ellipsoid behaviour is, in this case, formally equiva- 
lent to that of a biaxial rigid body, as far as the precessional 
motion is concerned. If 7o represents the moment of inertia 
of the spherical NS, the eigenvalues for the distorted NS are 
h= h = h- (l/3)AJf2 - (l/3)AJs and h = h + Alg, 
where AIn represents the centrifugal deformation. Since the 
latter term affects all eigenvalues in the same way it does not 
enter explicitly in the precession dynamics, which are deter- 
mined only by the magnetic deformation (see Appendix IX)) . 
Hence, the focus here and in the following will be on the 
magnetic deformation only. 

The spin energy of a rotating spheroid is minimized, 
at fixed angular momentum, when its moment of inertia is 
maximum. For a prolate figure, this is achieved when the 
symmetry axis is orthogonal to the spin axis. In the pres- 
ence of internal dissipative processes, the magnetic axis of a 
prolate NS will thus be driv en towards orthogonal rotation 
jMestel fc Takharlll97j : [jon cs 1976), which maximizes the 
efficiency of GW emissioi{3 (Cutleiil2002l ). 

For arbitrary GWs will be emitted at both the spin 
freq uency and its octave, w ith a total rate of energy emission 
(cfr. ICutler fc JonedboOll and references therein) 

E.W = -^^^ sin\{l + 15sin\) (1) 
c 

An orthogonal rotator emits GWs only at twice its spin fre- 
quency, at the rate given by eq. ([1} with x = 'i"/2 

E,„(x-V2) = -f ^^-^ (2) 

Given the generation mechanism for the superstrong inter- 
nal field, the magnetic axis is expected to be just slightly 
tilted, initially, to the spin axis and, for small x, the GW 
luminosity is largely suppressed (eq. [1} . Given the efficiency 
of magnetodipole radiation at birth, the orthogonalization 

^ A dipolc field induces an oblate deformation with the magnetic 
(symmetry) axis corresponding to the largest axis of inertia. In 
the presence of dissipation, the figure will thus be driven towards 
aligned rotation and no GW emission will result. 



process must thus be quick enough for strong GW emission 
to ensue promptly, in order to be competitive with magne- 
todipole radiation. 

Strictly speaking, the spindown luminosity of a mag- 
netic dipole rotating in vacuo depends only on the mag- 
nitude of the dipole component orthogonal to the spin 
axis. However, according to the standard pulsar model, 
an aligned rotator is expected to have a spindown lumi- 
nosity comparable to that of an orthogonal rotator. Stud- 
ies on the structure of force-free NS magnetospheres sug- 
gest that this is indeed the case, within a factor of order 
unity ( Contop oulos. Kazanas fc Fendtlll999l : [G"ruzinovir2006l : 
[Spitkovskv 2006). Based on these results, we assume a newly 
formed magnetar to have the magnetodipole luminosity of 
an orthogonal rotator. Further discussion of these issues is 
delayed to § [l] 

We neglect here the effects of strongly magnetised 
winds from newly formed magnetars spinning at ~ 1 ms 
period and restrict our attention to external dipole fields 
^ 5 X 10^'* G. Although such winds can be extremely ef- 
ficient in carrying away angular momentum (and spin en- 
ergy) from the NS in just a few minutes, their efficiency 
is expected to be negligible for dipole fields weaker than 
(6-^ 7) X 10^^ G (mompson et al.ll2004l : iBucciantini et all 
l2006l : lMetzger et al.ll2007f K 

To summarize, we describe the early spin evolution of a 
newly born magnetar as being driven by both magnetodipole 
and GW torques. Hence, as the tilt angle x increases, a 
sufficiently large (prolate) deformation and rapid initial ro- 
tation can make GW emission dominate angular momen- 
tum and rotational energy losses. For a given spindown 
torque, uj — — KcOj", the corresponding spindown timescale 
is Tad = uj/{2uj). For convenience, in the following we express 
the relative strength of the two torques at birth in terms of 
the ratio: 

where A = Kd/Kgw; GW losses are dominant for x > 1. 
The complete spin-down equation we adopt is thus: 

(Jj = ---^ Lo'-' - ^ LU^' = -Kd id^ - Kg„ LU^ (4) 

3 Ic^ Ilu 
where Eg„ is given by ([T}, or by ((2)1 when x ~ ^/2. 

2.1 Damping of freebody precession and 
orthogonalization 

We consider here in some detail the very first stages when 
the prolate NS is formed, its symmetry axis being just 
slightly tilted to the spin axis. We focus in particular on 
those processes that might affect the orientation of the sym- 
metry axis and thus promote (or prevent) prompt and strong 
GW emission. 

As stated above, freebody precession of the prolate 
spheroid is eventually damped by internal viscous torques, 
which redistribute angular momentum (with no loss) inside 
the NS so as to minim iz e the spin energy. Following the 
analysis of lCutleij (|2002l ). [stella et all (|2005| ) assumed that 
internal (viscous) dissipation of free p recession occurrs on a 
(very short) timescale Tart ^ Pp.eo l|Alpar fc Saulslll988l ). 
where Pprcc (< few seconds) is the free precession period. 
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This estimate results from the crust-core coupling caused 
by the interaction between superfluid neutron vortices and 
relativistic electrons in the NS core, where only electrons fol- 
low the instantaneous rotation of the crust. The application 
of this prescription to a newborn NS is subject to important 
caveats: neutron pairing in a sta te in t he core occurs 
at a ter nperature T,o„d < 2 x 10^ K (|TD96| and references 
therein, IPage et alj |2004 ) and crust formation also occurs 
at a tempe rature ~ a few x lO'^ K . The coupling mechanism 
studied bv lAlpar fc SauQ l|l988l ) does not apply when T > 
Tcond, since there is no superfluid and, likely, not even a 
proper crust. The NS is more like a self-gravitating, rapidly 
rotating, fully degenerate fluid mass. A proper account of its 
early cooling is thus of crucial importance in this context. 

If direct Urea processes occurr in the densest parts of 
the core, then cooling is extremely fast due to the very large 
neutrino luminosity. The NS temperature drops to 10^ K in a 
matter of minutes, as op posed to ~ 1 yr in the case of m od- 
ified Urea cooling (cfr. IPaee. Geppert fc Webeij l2006l and 
references therein). Therefore, if newly born magnetars cool 
through direct Urea processes, crust formation and tran- 
sition to superfluidity in the core occur ver y quickly. The 
crust- core coupling mechanism described bv lAlpar fc SauQ 
l|l988| ) could thus operate soon after NS birth and lead to 
the very short orthogonalization time estimated above. 

In the opposite limit, when NS cooling is driven by only 
the modified Urea reactions, the evolution is more complex. 
The temperature in t his case evolves as JOwen et al.l flOOSl 
and references therein, iPage. Geppert fc Webe rl l2006l ) 



T(t) 
109 K 



1_ f 10'^ K 

Tc \ Ti 



-1/6 



(5) 



where Tc ~ 1 yr and Ti ~ 10^" K is the initial tempjsrature 
of the NS at the end of the a — O, dynamo phase (! TD96h . 
As we are interested in the cooling at T > Tcond , no further 
neutrino-emitting reactions are considered, such as those oc- 
curring atjrwJTcondJn the so-called "minimal cooling sce- 
nario" (|Page et al.ll200j and references therein). According 
to eq. (O, the NS temperature will reach a value 2 x 10^ 
K in about 5 days after formation. During this time, the 
crust- core coupling mechanism introduced bv lAlpar fc Saulsl 
1 19881 ) cannot operate. In the absence of other viscous pro- 
cesses, the prolate spheroid will thus not be orthogonalized 
and GW emission will remain highly suppressed during those 
early days. Therefore magnetodipole radiation by a 10^* G 
external fleld will carry away most of the initial spin energy, 
thus spoiling the possibility of significant GW emission also 
at later times. 

However, the prolate NS is subject to at least two fur- 
ther (and competing) mechanisms that can in principle al- 
ter the orientation of its symmetry axis. First, GWs will 
be emitted even for a small initial tilt angle Xij though at 
a small rate. The corresponding radiation reaction torque 
will cause the spin and magnetic axes to align, regardless 
of wh ether the spheroid is oblate or prolate (|Cutler fc Joned 
120011 ). This process always acts so as to quench the GW 
emission efficiency. 

On the other hand, freebody precession induces internal mo- 
tions in the newly formed, fluid NS, that are required by the 
condition of hydrostatic equilibrium l|Mestel fc Takharl 19721 : 
|joneslll976l ). Dissipation of these internal motions through 



bulk viscosity is potentially able to orthgonalize the symme- 
try axis of the spheroid (the magnetic fleld axis) relative to 
the angular momentum vector, thus increasing the efficiency 
of GW emission. 

The relative strength of these two mechanisms will thus de- 
termine the early evolution of the angle x and the possibility 
of newly formed magnetars to become efficient GW emit- 
ters. Below we study these two processes in more detail, 
showing that orthogonalization through bulk viscous damp- 
ing is expected to always prevail on radiation reaction, in 
the parameter range of interest to our work. 

2.1.1 Equation of state and related quantities 

For our aims, the equation of state (EOS throughout) of 
NS matter enters through its role in determining the global 
structural properties of the NS, such as the mass-radius re- 
lation and the moment of inertia. We dot not consider the 
detailed microphysics that determine the EOS; rather we 
adopt a phenomenological approach, parametrizing all re- 
sults in terms of the NS mass and r adius. 
As shown by iLattimer fc PrakashI (1200 ll ) most NS EOS's 
are well approximated by a polytrope of index n — 1, 
with P = kp^, as far as their global properties are con- 
cerned. Here k is determined by the stellar radius (see Ap- 
pendix |B]|. We also adopt the approximate formula for th e 
NS moment of inertia given bv ILattimer fc PrakashI l|200lf ): 
O.SSMR^ ~ 1.4 X 10*''(M/1.4 M(7) )(R/12 km)^ g cm^ 
Finally we note that, as discussed in ILattimer fc PrakashI 
(2007), internal magnetic fields weaker than ~ 10^* G are 
not expected to sizeably affect the microphysics that deter- 
mine the EOS. 

For the sake of completeness, we also report here the 
resulting relation between the measured P and P of a NS 
and the corresponding magnetic field strength as derived by 
the magnetodipole spindown formula (see eq. |3| . 
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2.1.2 The Magnetically-Induced Distortion of a 
Neutron Star 

One of the key parameters of interest is the deviation from 
spherical symmetry (es) of the NS, that is caused by the 
anisotropic pressure of the internal magnetic field. This de- 
termines the frequency of freebody precession and the GW 
luminosity of the rapidly spinning NS. In general, es will 
be determined by the volume-integrated ratio of the mag- 
netic to gravitational binding energy densities, times some 
numerical coefficient accounting for field geometry and NS 
structure. Setting the latter factor to unity, and assuming 
an approximately constant density in the NS core 
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as deter mined from g eneral arguments based on the virial 
theorem (|Cutlerll2002l and references therein). In the above 
definition, the integral is extended to the whole NS volume 
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where the magnetic field is present. We are assuming the 
internal field to be mostly toroidal, the poloidal component 
being sufficiently small to be energetically negligible, as en- 
visaged in the millisecond magnetar formation scenario. Eg 
in eq. (O is the gravitational binding energy of the NS that, 
for a polytrope of index n = 1, is Eg = (3/4)GMVi?. 
The value of es can thus be approximated as 



es ^ -1.15 X 10-^ f -^^) ( (- 



M 



AUr. 



(8) 



that is negative for a prolate shape. The volume-averaged 
strength of the internal (toroidal) magnetic field is: 



Bt 



2 X IQis G 



93 ( y/y R \ 

VlO«' crg^ Vl2 km^ 



-3/2 



(9) 



Uncertainties in the exact value of e_g as a function of the 
(unknown) internal field distribution can be parametrized 
by adding a further multiplicative facto]0 (rj) in the right- 
hand side of eq. ([7]) where, in general, 77 > 1 is expected. We 
neglect this factor for clarity, but stress that the quantity 
Eb used throughout can more generally be substituted with 
rjEs- We explicitly assume rj = 1 from here on. 



2.1.3 Gravitational radiation reaction 

The effect of gravitational radiation reacti on on damping 
of free body precession was investigated by ICutler fc Joned 
who showed that GW emission always drives the tilt 
angle x to zero, independent of whether the NS is oblate or 
prolate. In the limit of small x, these authors express the 
freebody precession damping timescale {t„) as: 



5c" 



d(sin x)/dt 
E 



3.7 I 



4 / R N -4 



VlO^Oerg/ Vms/ Vl2km/ 



(10) 



The importance of this effect on the early evolution of a 
newly formed magnetar can be determined by comparing 
this timescale to the bulk viscosity dissipation timescale and 
to the magnetodipole spindown timescale. 



^ The exact value of eg is uncertain as it depends on the mag- 
netic field configuration and EOS. Using a general rela t ivistic 
treatment of the NS structure, iBonazzola &: Gourgoulhoni (Il996l ) 
studied the way in which the numerical factor in eq. jTj changes 
for different configurations of the internal magnetic field. These 
authors showed that this factor can be higher than 15/4, or even 
much higher, for specific magnetic field geometries (for example 
the case of an internal field that is mainly concentrated in an outer 
shell of the NS core, or the case of an intermittent fi eld distribu- 
tion rather than a uniform one) . iHaskell et al.l 1I2OO8I) studied the 
same problem allowing for different EOS's and combinations of 
the toroidal and poloidal components of the internal field. Their 
results confirm that the magnetically-induced ellipticity can be 
in general somewhat larger, or significantly larger in particular 
cases, than the estimate given by eq. ([7]l and 

Overall then, for a given magnetic energy, larger ellipticities than 
given by eq. (0 are a possibility worth of further investigations. 



2.1.J^ Bulk viscosity in newly born NSs 

As shown bv lMestel fc Takhaij l|l972l) . a field of internal mo- 
tions is excited in order to maintain hydrostatic equilibrium 
everywhere within a fiuid star undergoing freebody preces- 
sion. These motions are periodic, with the same frequency as 
freebody precession and amplitude proportional to the (non- 
spherical) centrifugal deformation of the star (pn oc n^) and 
tilt angle x of the magnetic axis. In this section we discuss 
the damping of such motions by bulk viscosity and the asso- 
ciated growth of the tilt angle x, within a non-superfiuid NS 
made only of neutrons, protons and electrons {npe matter). 

The assumption of pure npe matter corresponds to the 
least favourable case for damping through bulk viscosity and 
likely not the most realistic. At the highest densities of NS 
cores muons and possibly hyperons (A, E~) are expected 
to appear, significantly increasing the bulk visc osity coef- 
ficient fe.g. |Jones|[l97g : iLindblom fc Qwenll2002l and refer- 
ences therein). In this respect, our calculations are conser- 
vative. Detailed models should lead to significantly shorter 
damping times for the excited modes. 

Fluid bulk viscosity is a consequence of the finite time it 
takes for the fluid to react to changes in one of its thermo- 
dynamic parameters, adjusting the others to their new equi- 
librium values. In hot NS matter, displacement of a fluid 
parcel from the equilibrium position causes departure from 
chemical equilibrium between particle species: bulk viscosity 
is thus determined by the activation of /3-reactions trying to 
restore chemical equilibrium. A general expression for the 
bu lk viscosity coefficient in N S matter (^) has been derived 
bv lLindblom fc Owed l|2002t ) 

nr {dp/dx)^dx / dn n{dp/dx)^dx/dn 



Re(C) = 



(11) 



1 -I- {wry Lj^r 

where n is the particle density, p is the pressure (whose 
derivative is calculated at constant proton fraction, x = 
np/n„ ~ Up/n) and x is the equilibrium value of x. 
The quantities in the denominator are the mode frequency 
(ojpro ~ EBHspin) and /3-reaction timescale (rp). Defining 
the characteristic damping time of the e xcited motions as 
Td = 2Eprc/|Epro| one can express it as fcf. lOwen et al.iri998l l 



Td = 



/ C|V • 5v\^ dV 



2 Ep 



2 Ep 



.(12) 



where 5v is the velocity perturbation associated to the mode. 
In the last step we have u sed the relation V -Sv = iuj{An/n) 
(|Lindblom fc Owen|[20o3 ) and substituted the Eulerian per- 
turbation, Sp, to the Lagrangian one, Ap. A word of cau - 
tion is needed here. The analysis bv lMestel fc Takharl (|l972l ) 
gives V ■ ^ = 0, i.e Ap — and Sp = — ^ • Vp, a condition 
that would imply the absence of bulk viscosity. However, this 
result is a direct consequence of having considered strictly 
adiabatic fluid motiouijj within a radative (sub-adiabatic) 
layer, in a first-order perturbative analysis. As opposed to 
this, a polytropic NS would have an almost adiabatic gra- 
dient, in which case the calculations by Mestel fc Takhar 
would leave V • ^ totally unconstrained, revealing the need 
for a higher-order pertubation analysis. Indeed these authors 

In fact, bulk viscosity is a result of deviations from strict adia- 
baticity of perturbations. 
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discuss the (non)applicability of their conclusion to a con- 
vective (almost adiabatic) zone, at the end of their section 3. 
They discuss a number of approximations in their analysis 
that may well break in the case of a very rapid rotator (such 
as a millisecond magnetar) thus introducing non-negligible 
higher-order terms in the perturbations. In these cases, the 
internal dynamics of the oblique rotator becomes much more 
complicated and a detailed analysis is beyond our scope here. 

In general, however, motions of the same order o f mag - 
nitude as the one discussed by iMestel fc Takharl (|l972 ) 
are always expected to occur, whose rate of expansion 
(V ■ must be calculat e d exp licitly. Here, in analogy to 
iReisenegger fc GoldreichI (|l992 ) and to several treatments 
of bulk viscosity for r- modes i Lindblom. Owen fc MorsinkI 
1 19981 : lOwen et al.lll998l ) , we assume that Ap would be of the 
same order of magnitude of 5p, leaving more detailed analy- 
ses to fut ure wo rk. To stress the importa nce of this aspect we 
note that ILindblom. Mendell fc OwenI l|l999l ) carried out a 
second-order analysis for r-mode damping through bulk vis- 
cosity and derived an order of magnitude longer dissipation 
timescale than those based on the assumption Ap ~ 5p. 

We finally note that Sv ~ eBOjR ~ 10*^ cm s"'^, while the 
Alfven velocity is va = B j ^f^Kp ~ (10*^-^10^) cm s~^ for the 
parameters of interest. Alfven waves are thus very efficiently 
in maintining rigid rotation of the fiuid star, despite the 
perturbation in principle introduced by ^. 

The calculation of the perturbation amplitude {5p/p) 
is detailed and discussed in Appendix |B] We report here 
our result for the damping timescale of freebody precession 
according to eq. (fT^ 



Td 



13.5 



cotan X / Eb 
1 + 3cos2x V lOSOerg 



\ ms / 



M 



1.4 M 



T 



IQiOK 



(13) 



a fairly short time even for small initial values of the tilt 
angle (x»). 

Given ra , we can eventually calculate the growth time of 
the tilt angle, r^. As shown in Appendix|X]damping of free- 
body precession reduces the NS spin energy, at a constant 
angular momentum, by changing the tilt angle x. From eq. 
we get: 



Td = 



cosx 



2 

COS x 



sinx 



X sinx sin'^x d(sin x)/dt 



— cotan^x T^ 



(14) 



Here, is the timescale to be compared to (eq. llOp in 
order to determine the relative importance of bulk viscosity 
and radiation reaction in the evolution of x- In the limit of 
small tilt angle, for which cosx — 1, we have 

zi ^ 10-5 ( \' fi_\ f^y' ns) 

VlflSOergy Vms^ KIO^^'k) \l2Km) ^ ' 

From this we derive the condition for bulk viscosity to 
largely prevail on gravitational radiation reaction, so that 
orthogonalization is essentially unaffected by radiation- 
reaction. For definiteness we require < O.It^, a condition 
that translates to 



V10"'K/ Vms/ V 



2/3 /12km\4/3 



R / 



erg 



(16) 



oiQb < 9 X W (T/10"'K) (P.^s)'''^ G. 
Soon after formation, the temperature T decreases on a very 
short timescale, not much different from itself, even for 
the slow cooling given by eq. The angle x a-nd T thus 
evolve on comparable timescales and just considering the 
damping timescale at a given temperature is not appropri- 
ate. Rather, the coupled evolution of x and T must be solved 
in order to derive a reliable estimate of the time it takes for 
X to grow to large values. 

Inserting the cooling history ([5]) in (|13|l or (|14|l . the 
resulting equation for x can be solved with initial con- 
ditions X = Xi and Ti = 10^° K. Since Td(t) = 
N T~^{t) cos2x/[sin2x(l -I- 3cos2x)], the expression T'f (t) 



[(t/30) + 1] gives: 



cotanx 
1 -I- 3cos2x 



dx- 



dt 



N 



(17) 



whose solution is: 



sin^x 



1 -I- 3cos2x 



sin^Xt 
1 4- 3cos2xi 

sin^X» 
1 -I- 3cos2xi 



f- 

V30 



240/N 



-1440/N 



(18) 



From the above we can obtain the time (in seconds) or the 
temperature (in units of 10^" K) at which a sufficiently large 
value of the angle x is reached starting from a given, small 
initial tilt angle Xi- 

The requirement that the damping time of free precession be 
much shorter than the initial spindown timescale of the NS 
allows two key constraints to be met jointly and our scenario 
to maintain full self-consistence; first, damping of free pre- 
cession will be well described in terms of the approximation 
of constant angular momentum. Second, efficient GW emis- 
sion will ensue quick enough for the NS initial spin energy 
to be still fully available. 

The initial timescale for electromagnetic spindown is 
rj„ ~ 1.1 Pi,B^i4 d, and in a time ~ Clrj^ less than 10% 
of the initial spin energy is lost to magnetic dipole radia- 
tion. Hence, we consider this as the longest time over which 
orthogonalization must take place for our scenario to apply. 
Therefore, given x* = 1 deg, the angle x will grow to a suf- 
ficiently large value - say, x = 60 degj - in a time shorter 
than 0.1 ri„ i£| 
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Bd,14 
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Ml.. 



Bd.l4 R12 



(19) 



Note that the numerical coefficient in the last step is ~ 5 if 
Xi = 2 deg, and ~ 2.8 if Xi = 0.1 deg. 



For small x the temperature is always very near to 10"'^'' K. 
Therefore, this constraint is easily met if B is not very close to 
IQl^ G. 

5 For X — 60 deg GW emission is very efficient and can be ap- 
proximated by eq. ((2]l 

^ Note that for any Bd > 6x lO^^ q ^nd Pi,,„, ^ 3 ms, 0.lTi_em ^ 
3 days. The core temperature will accordingly be > 2 X 10^ K, 
so that our condition covers essentially the range of temperature 
before the core becomes superfiuid. 
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Figure 1. The region of parameter space, in the external dipole 
field strength (Bd) vs. internal magnetic energy (Eb), where the 
timescale for the angle x to grow to a large value (60 deg) is less 
than one tenth of the initial spindown timescale through magne- 
todipole radiation. The favourable region lies on the left of the 
corresponding limiting curve. For each value of the initial spin 
period, two curves are plotted for two different values of the ini- 
tial tilt angle, Xi = 1 (left) and 2 deg (right). Note that further 
decreasing Xi from 1 deg to 0.1 deg shifts the curves to the left 
by a slightly larger amount than the decrease of Xi from 2 deg 
to 1 deg (cfr. eq. I19I I. This would still leave ample room in pa- 
rameter space for the fastest spinning magnetars, while it would 
drastically reduce the available parameter space for magnetars 
with Pi > 2.5 ms. 



Eq. ([19} is represented in Fig. [T] by two curves corre- 
sponding to initial spin periods of 0.97 and 2.58 ms, that 
approximately bracket the relevant range of spin periods for 
newly formed magnetars. Each curve is plotted for two dif- 
ferent values of the initial tilt angle, Xi = 1 s^nd 2 deg (see 
caption for further details). 

The constraint shown in Fig. [T] will be used, together 
with independent constraints derived in later sections, to 
identify the region in parameter space where GW emission 
efficiency from newly formed magnetars is optimized. We 
stress again here that all our calculations were carried out 
by assuming pure npe matter, a conservative assumption 
that minimises the efficiency of bulk viscosity. Finally, in 
Appendix |B] we calculate the centrifugal distortion of the 
density profile in the slow-rotation limit, which is also likely 
to underestimate the deformation of a millisecond magnetar. 
This translates into an underestimate of 5p in eq. (|12|l and, 
thus, of the efficiency of bulk viscosity. 

We conclude that, despite having chosen a "worst case 
approximation", damping of freebody precession through 
bulk viscosity in a newly formed, rapidly spinning and 
strongly magnetized NS can be very efficient in the param- 
eter range considered here. Therefore, strong GW emission 
from an almost orthogonal, rapidly rotating NS is likely to 
ensue quickly as a consequence of a strong toroidal magnetic 
field and of the efficient dissipation of its freebody preces- 
sion energy. The implications of this are explored further in 
the next section. 



AMPLITUDE AND DETECTABILITY OF 
THE EMITTED SIGNAL 



IStella et all (|2005l ) calculated the expected signal-to-noise 
ratio (S/N) that a putative GW signal from a newly born 
magnetar would have, for an optimal (matched-filter) detec- 
tion, adopting the broadband design sensitivity of Advanced 
Ligo. At frequencies between 0.5 and 2 kHz this is well ap- 
proximated by 



S;.(/)^So f 



(20) 



where Sh{f) represents the one-sided spectral noise distri- 
bution of the detector, / = 2a; the GW signal frequency 
and So ~ 2.1 x lO"''^ Hz"^ (|Owen fc Lindblom|[20o3 : ICutlej 
l2002i and references therein). Note that the designed sensi- 
tivity for A dvanced Virgo is very similar, in this range of 
frequencies (jLosurdol |2007| ). so that our calculations hold 
essentially for both detectors, frequency, f — uj/n for an or- 
thogonal prolate rotator. 

We re-address here this point to better qualify the role of the 
NS ellipticity in the detectability of the signal. We also cor- 
rect a (sma ll) numerica l error in the calculated S/N curves 
in Fig. 1 of IStella et all (|2005l ), whose conclusions maintain 
their general validity. The S/N of an optimal signal search 
is defined as 



S/N = 2 



Shif) 



1/2 



(21) 



Here h{f) is the Fourier transform of the instantaneous sig- 
nal strain h[f{t)] that, i n the stationary phase a pproxima- 
tion, is expressed as (cfr. lOwen fc Lindblomir2002l and refer- 
ences therein) 



(22) 



where the time derivative of / is obtained from eq.(|4]), 
since / = ui/tt. We adopt the expression for the strain 
amplitude - averaged over source or ientation - given bjQ 
lUshomirskv. Cutler fc BildstenI (|2000l ): 



1/2 



Gles 
Dc4 



(23) 



where D is the source distance. Further averaging over the 
detector antenna pattern, eq. (|2ip gives the optimal S/N 
ratio: 



S/N = 



Sh{f) d//dt 
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(24) 



where we have set a — Kd/(7r^Kg„) — Aj-K^ . Substituing the 
numerical values: 



6 (^j(,M)erg) (ti^) (t5hi) (rm^) 



-1/2 



i D V 
\20Mpc j 



21nf +ln4^ 



1/2 



(25) 



^ Comparin g with the "optim ally-oriented" strain amplitude 
fto given by lAbbot et aO l|2007l V we obtain the relation ha = 
4/(5v^)ho. 
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Figure 2. The optimal (matched-filtered) signal-to-noise for searches of the GW signal from newly formed magnetars in the Virgo 
Cluster with Advanced Ligo /Virgo, calculated through eq. (|25| l. Four different values of the initial spin (indicated in the figures) are 
considered and results are shown as contour levels (at the values of S/N indicated on each line) in the Ba vs. Eb plane. 



In Fig. [2] we show the curves of constant S/N in the B d vs. 
Eb for selected values of the initial spin, as derived from eq. 
I|25p . More details are given in the caption. According to eq. 
I|24p the maximum S/N is obtained in the limit a — > 0, which 
depends only on the initial spin energy of the NS and not 
on its oblateness. Clearly, this value [S/N] max is attained as 
the magnetodipole spin-down torque disappears, so that all 
of the initial spin energy of the NS is lost to GWs. This can 
also be seen directly by substituting the expression for the 
pure GW-driven spindown in eq. (|24|l . which gives 



[S/N], 



4.5 



D \- 



(2OMPC, 
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12Kri 
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(26) 



3.1 Quantifying the gravitational and 
electromagnetic energy output 

Thus far we have discussed the conditions under which the 
hypothesis that newly born magnetars be detectable sources 
of GWs with next generation detectors can be true. X-ray 
observations of Supernova Remnants (SNRs) around mag- 
netar candidates in the Galaxy give us clues on the actual 
viability of this hypothesis (cf. §[l|. In this section, we show 
that our scenario is indeed fairly consistent with such ob- 



servations. We find that, if magnetars at birth were de- 
tectable sources of GWs from Virgo-cluster distance, then 
SNRs around them would likely show no significant excess 
of energy injection with respect to "ordinary" SNR^ 

We begin calculating the total integrated energy emit- 
ted via GWs by a newly born magnetar as: 



dt = 



du) 



(27) 



For the spindown model in eq. (Q, the above integral gives 



Jo 



= Egpin.i 



duo = I 
1 - 



(28) 



Eq. p8|) expresses AEg„ as a fraction {5) of the initial 
spin energy of the NS. The remaining energy AEc, = 
Espin,i(l — 5) is the maximum that can be transferred to 
the SN ejecta. Next, we define the quantity AEtr = /3AEo„ 
as the amount of that energy that is effectively transferred 
to the ejecta. In the absence of additional competing torques 
and/or a physical prescription for the transfer mechanism, 
we will assume it to be perfectly efficient (/3 = 1), although 



" Assuming GWs do not appreciably transfer energy to the Su- 
pernova ejecta 
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leaving the explicit dependence on P in the calculations that 
follow. 

Recalling the definition of x as the ratio of gravitational to 
magnetic dipole torque at birth (eq. [Sjl and using eq. (|28p 
we then get 



ln(l + x) _ AE„ 



1 AEt, ^ 1 

/3 Eapin^i /3 Ejipi, 



SNR 



(29) 



where iSg^^ is the amount of energy tra nsferred to the 
ejecta as determined from observations (jVink fc Kuipeij 
I2OO6I ). Observations provide an upper bound to this quan- 
tity, from which the inequality in the last step in eq. 
(|29p obtains. We derive constraints on magnetar param- 
eters at birth by usi ng the upper bound determined by 
IVink fc Kuip"^ i 20061 ) discussed in this paper and our GW 



emission scenario. 



First fix a value for x: eq. (|29p gives the maximum value 
of the initial spin energy or, equivalently, )„,ax compatible 
with the upper limit on Eg'JiJj^. Accordingly, the maximum 
value of A = {cOi)i„ax/x is constrained (cfr. eq. |3}: 



(31 ln(l +x) 



(30) 



Inserting eq. (|30|) in the definition of A= Kd/ Kg„ (cfr. § 
[2| we get the following condition, for the electromagnetic 
output of a newly born, millisecond spinning magnetar to 
be always ^ 10^^ erg: 



G l3Gln(l+x) 3/2 
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(31) 



Therefore, given x one obtains uji and, eventually, the ap- 
propriate curve in the Bd vs. Eb plane. Results are shown in 
Fig. [3] where curves in the Bd vs.Eb plane have been drawn 
for four different values of x = 10, 22, 100 and 150, corre- 
sponding to the initial spins indicated in each plot. Curves 
of given signal-to-noise ratio from Fig. ([2]) are also plotted 
together with curves defining the region where the orthogo- 
nalization timescale is sufficiently fast (see Fig. [ij . 

We conclude that there exists a wide region of the pa- 
rameter space where all constraints for efficient GW emis- 
sion from newly born, millisecond spinning magnetars are 
met jointly. 



4 ELECTROMAGNETIC SPIN-DOWN AT 
BIRTH 

The ideal magnetodipole spindown formula (cfr. eq. H} im- 
plies that a sizeable fraction of the initial spin energy of a 
millisecond spinning magnetar would be lost to electromag- 
netic emission, for Bd ^ 3 x lO'^* G. 

SGRs, and two out of seven AXPs, have spin periods 
and period derivatives leading to esti mated dipole fields of 
a few times lO" G, up to ~ lO^'^ (cfr. IWoods fc ThompsonI 
I2OO6I and references therein), some of which are apparently 
out of the optimal range for GW emission at birth, if they 



maintained their dipole field since then. Although relation 
((4]) is certainly true to a good degree of approximation, the 
ideal magnetodipole formula does not appear to hold exactly 
in the few isolated NSs for which sufficiently accurate timing 
measurements exist. 

Observations show that the spin-down of pulsars, as- 
sumed to be wholly due to an electromagnetic torque, might 
be a somewhat weaker function of the spin frequency than 
expected in the ideal case. Indeed, in the few pulsars where 
the second time derivative of the spin frequency (i>) could 
be measured, a "braking index" n = (uiy/u^) < 3 is typi- 
cally deri ved. In particular, the Crab puls ar has a measured 
n = 2.5 (iLvne. Pritchard fc Smithlll988l ) and PSR J1119- 
6127 has n = 2.9 |Camilo et al.1 l200d 'l. TLivingstone et al.1 
have recently measured braking indices of three rela- 
tively young X-ray pulsars. Again, all results give n < 3 (2.14 
for PSR B0540-69, 2.84 for PSR B1509-58 and 2.65 for PSR 
J1846-0258). The Vela pulsar has the lar gest known devia - 
tion from the ideal case, with n = 1.4 (|Lvne et al.lll996l ). 
although subtraction of the frequent glitches from its spin 
history presents difficulties and that result should be taken 
with some caution. 

Similar (or even greater) unc ertainties are associat ed to the 
putative n ~ 6 determined bv lMarshall et al.1 l|2004h for PSR 
J0537-6910 (the "Big Glitcher"). This source is found to 
exhibit very frequent glitches; post-glitch "recoveries" are 
known to significantly affect the second time derivative of 
the spin frequency. For this reason, the measured values 
of i> from such a frequent glitcher are probably unsuitable 
to determine reliably th e "true" secular spindown trend. 
iMiddleditch et al] (I2OO6I ) discuss in detail several aspects 
concerning the timing properties of this source, arguing for 
a much smaller value of n. 

The examples of 10'^ -r- 10* yrs old pulsars with brak- 
ing indices lower than 3 suggest considering the same pos- 
sibility for magnetar candidates, whi ch indeed have similar 
ages and similarly high glitch activ ity l|DairOsso et al.ll2003l : 
llsrael et all l2007l : iDib et al.1 [20081 ). No measurements of n 
has yet been obtained for AXPs and SGRs. It would be in- 
teresting to know from the data whether these sources share 
this same spindown property of other isolated NSs. 

Several mechanisms able to produce a braking in- 
dex smaller than 3 have been proposed in the literature, 
since the early measureme nts of n in the Crab pulsar. 
iBlandford &: Romanil (|l988i ) discuss in general the effects 
of evolving the parameters of the magnetodipole torque. 
They identify a secular increase in the surface magnetic 
field, at least over the first ~ 10"^ -r- 10* years, as a 
most promising explanation for n < 3 in young NSs. 
Amplification of crustal fields through a thermomagnetic 
battery over this early phase had been previously dis- 
cussed (|Blandford. Applegate fc Hernauistlll983l and refer- 
ences therein), which could find application also in this con- 
text. 

Alternatively, several suggestions have been made for 
a secular increase of the tilt angle between the mag- 
netic dipole and the spin a xis, based on a wide vari- 
ety of physical mechanisms |Goldreicbl Il970l: iRuderman 



199g 



199ll : lLink. Epstein fc Bavmlll992l : lRuderman. Zhu fc Chen 
More recently, the possibility that the zone of closed 



field lines in the magnetosphere (the corotating region) does 
not reach the light cylinder has been widely discussed. 
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Figure 3. Curves defining the region of parameter space - in the Bd vs. Eb plane - where all constraints for efficient emission and 
detection of GWs are fulfilled, for four initial spin periods. The region on the left of (and below) the solid curve labelled r^^ < 0.1t^'„ 
is where the tilt angle x becomes 60 deg in less than one fifth of the initial magnetodipole spindown time. The region on the right 
of (and below) the straight, solid curve is where the energy emitted electromagnetically by the newly formed magnetar is ^ 10*^ erg. 
It is obtained from eq. I I31I I with four different values of x = 150, 100, 22, 10 from top to bottom and right to left, corresponding to the 
initial spins indicated in the figures. The region on the right of the dashed (or dot-dashed) curves is where a match-filtered search with 
Advanced LIGO/Virgo would give a signal-to-noise ratio greater than indicated in the figures, with the source at the distance of the 
Virgo Cluster (20 Mpc). 



both in t he general context of stud ies of pulsar magne- 
tospheres (IContopoulos. Kazanas fc Fon dt 1999j Gru zinovi 
l2006l : IContopoulos fc Spitkovskvl bood : ISpitkovskvl I2OO6I ) 
and in relation to the magnetosphe r ic str u cture of newly 
born magnet ars (jBucciantini et al] I2OO6I : iMetzger et al.l 
l200'<t) . 

Simulations by ISpitkovskvl (|2006l ) have shown that the 
closed field line region reaches the light cylinder in a matter 
of one (at most) rotation period, as a model NS is set into 
rotation. However, it is not obvious whether the region of 
closed field-lines can subsequently track the expansion of the 
light cylinder as the NS spins down or, rather, lag behind it 
at an increasing distance, which wou ld naturally lead to n < 
3 IContopoulos &: Spitkovskvll2006l and references therein) . 

Without indicating any particular mechanism, we dis- 
cuss the consequences that an electromagnetic braking index 
n < 3 would have, at a purely phenomenological level. 
Considering a generic model for electromagnetic spindown 
with braking index n, one has 



■,(") 



(32) 



where the constant K^"^ = (tJo/'^o) ^nd the subscript "0" 
refers to present-day values of the parameters. According 
to eq. p2|) . given measured timing parameters (ojctjo) of a 
source the expression for its spindown at birth in the case 
n 7^ 3 and case n = 3 are related through 

(.) ^ ^ _j,,3)^3 ^ ,3) (^y^ (33) 

where the subscript "i" indicates quantities at birth. 
Note that the quantity in parentheses is smaller than unity 
for n < 3. It is therefore natural to ask what value of n 
would be required in magnetar candidates (AXPs/SGRs), 
given their measured loo and luq, for their spindown at birth 
to have been dominated by GW emission rather than by the 
magnetic dipole radiation. 

In order to answer this, we define the maximum allowed 
strength of the electromagnetic spindown at birth (tji,max). 
Given the results of the previous section, this maximum 
value will equal that the ideal magnetodipole formula (n =3) 
would give for Bd = Bd(max) = 2x 10^^ G (a value that guar- 
antees strong GW emission at birth, cfr. § 13. ip . Therefore 
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we write, uji 



K 



(3) 



of and the last step of eq. (|33p 



must be smaller than [i;i,max, which gives: 
Kd — s? ' (max) 

From this, since the ratio of the torque functions corresponds 
to the ratio of the magnetic dipole fields: 



(34) 



(3 - n)Log ( — 

\LJO 



> 2 Log 



^ 3-2 



Bd 



Bd (max) 
Log [Bd/Bd(max)] 
3 + LogPo — LogPi,, 



(35) 



Among SGRs and AXPs, SGR 1806-20 has the strongest 
inferred dipole f ield (~ 1.1 x lO^'^ G with R= 12 km and 
M= 1.4M0, cfr. IWoods fc Thompsonll2006h . thus requiring 
the largest deviation from n = 3 in the hypothesis discussed 
here. Even assuming a (relatively) slow initial spin period 
for this source. Pi = 3 ms, eq. (|35p gives n ^ 2.6, a wholly 
plausible value compared to other isolated NSs. Note that 
the constraint is slightly weaker for other SGRs and/or con- 
sidering a spin period at birth shorter than 3 ms. For AXPs, 
the limit on n ranges from 2.7 to 2.85, for an (unfavourable) 
initial spin of 3 ms. This speculative argument would require 
direct measurements of braking indices in magnetar candi- 
dates. However, our aim here was to emphasize the depen- 
dence (strong, in some cases) of the calculations of previous 
sections on a number of poorly constrained physical param- 
eteres, and the importance of further studies on all of the 
above aspects. 



5 THE DECAY OF CORE FIELDS IN THE 10^" 
G RANGE 

The secular evolution (and dissipation) of the magnetic field 
in NS c ores wa s studied in detail by Goldreich & Reisenegger 
(1 992) (|GR92| from here on) and their analysis was extended 
bv lTD96l to the specific case of magnetar fields in the 10 G 
range. In magnetars, large-scale field instabilities lea ding to 
fast dissipation events were studied in detail as well ||TD95I : 
ITDOII : iLvutikovl |2003| ) . in order to interpret the powerful 
bursts and flares of SGRs. 

In GR92 three separate processes for secular field evolu- 
tion were identified two of which, ohmic dissipation and am- 
bipolar diffusion, are dissipative while the third. Hall drift, 
conserves magnetic energy. In particular, ambipolar diffu- 
sion was found to be more sensitive to the field intensity 
(GR92) which in fact implies that, while this process is not 
very important in normal NSs, it is the main mechanism for 
direct field decay in magnetars (TD96). 
Hall drift can affect indirectly the evolution and dissipa- 
tion of magnetic fields in NS interiors, on longer ti mescales 
than those char acteristic of ambipolar diffusion (|TD96| : 
I Arras et al.ll200^). As suggested in GR, 92 (and recently stud- 
ied in detail bv lCumming et aLll2004l ). excited Hall modes 
of field diffusion could decay (or cascade) to shorter wave- 
lenghts, that are subject to enhanced ohmic dissipation. 
This has particular relevance for accelerating fiel d decay in a 
magnetar's crust, as recent s tudies pointed out (|Pons et al.l 
120071 : H ons fc Geppert|[2007l ). Furthermore, Hall diffusion 
can drive - yet conserving the total energy - an initially 



stable MHD configuration close to a new equilibrium con- 
figuration with smaller total energy. A point can be reached 
where the field suddenly relaxes to the new equilibrium, if 
(sufficiently fast) fiuid motions are allowed within the stably 
stratified NS interior. 

As long as the early (ages much less than ~ lO** yr) evo- 
lution of magnetars is concerned, however, ambipolar diffu- 
sion in the NS core is expected to be the dominant mode of 
field decay. 

Ambipolar diffusion drives a slow motion of charged parti- 
cles with respect to background neutrons, which is opposed 
by both particle friction and chemic al pote ntial gradients 
in the stably stratified NS medium. IGR92| identified two 
separate modes of ambipolar diffusion, differing by their 
effect on chemical composition. The solenoidal mode does 
not perturb chemical equilibrium and thus is counteracted 
only by particle friction. The irrotational mode, on the other 
hand, does perturb chemical equilibrium and cannot evolve 
on timescales shorte r than the /3- reaction timescale. 
As shown in ItD96| . /3-reactions are very efficient at eras- 
ing chemical equilibrium imbalance when T > Tt, ~ 
5.73 X 10** (pi5/0.7)^/^^ K. At these high temperatures, 
both modes of ambipolar diffusion are effectively opposed 
by neutron-proton friction only. Field decay occurs on the 
same timescale in both modes (GR92) 



(early) 



AB2 



2.2 X 10' 



where L and a are the characteristic scale of variation of the 
Lorentz force and chemical potential, respectively, (LI a) ~ 
9.16 Tip~^^'-^{L/2 km) (GR92). The latter is > 1 at high 
temperatures, given the high efficiency of /3-reactions, while 
it becomes < 1 as the temperature drops and the efficiency 
of /3-reactions decreases. At T> 10^ K, the timescale ([36} 
is much longer than the NS age or its cooling timescale. 
Therefore, field decay is negligible as long as the temperature 
is this high. 

As NS cooling proceeds, the point is reached (at T ^ 
Tt,) where chemical equilibrium imbalance becomes the 
main obstacle against which magnetic stresses must work 
to drive particle diffusion. This affects only the irrotational 
mode, while the solenoidal mode still decays on the timescale 
of eq. (I36|l . Hence, the two modes grow at different rates 
with the irrotational mode evolving on a longer timescal^EI 
(TD96, GR92) 



{late) 
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1016G 



yr (37) 



TD96 considered in detail this lower-temperature regime. 
Based on stability arguments, these authors suggested that 
the solenoidal mode is expected to carry just a small frac- 
tion of the magnetic energy, most of it being tapped by the 



^ The two timescales are formally equal at T = Tt,, but the 
irrotational mode becomes quickly much slower below T = Ttr. 
For example, a 20% decrease of T below T = Tt, gives a 6 times 
longer decay timescale for the irrotational mode. 
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Figure 4. Expected evolution of the magnetic field intensity - at 
a specific value of the density p = 7 X fO^'* g cm~^ = 2.5pnuc - 
for three selected values of the initial intensity (Bi). Temperature 
equilibrium according to eq. 1 144 I I is assumed at each time. The 
tracks converge, as the field strength ap proache s 7 X fO^^ G, 
towards the asymptotic solution given bv lTD96i 



irrotational mode. The main conclusion of this scenario is 
that only a tiny fraction of the magnetic energy reservoir in 
the NS core is lost either in the high-T regime, or via field 
decay through the solenoidal mode at lower temperatures. 
Most of the magnetic energy dissipation occurs via the slow 
decay of the irrotational mode, at T ^ Tt,. This is based on 
the assumption that field decay, and the irrotational mode 
in particular, is effectively frozen at T > Tt,. 

The above scenario holds for fields in the 10^^ G range. 
According to eq. (ea. l36p . field decay at higher temperatures 
is significantly faster for stronger magnetic fields. Further, 
dissipation of even a small fraction of the magnetic energy 
reservoir may in principle affect NS cooling, if the reservoir 
is sufficiently large. In fact, field decay at T> 10^ K is not 
frozen if B is larger than 10^® G. In analogy to the treatment 
of TD96, we check here whether an equilibrium condition 
between heating and cooling in the high-T regime can apply 
as well, with field decay described by eq. (|36|) . For uniformity 
with that work we adopt the same normalizations for the 
parameters used bv lTDQq . 

We can write the heating rate per unit volume through 
field decay 



[earl] 



3.69 X 10' 



erg cm s 



(38) 



while the cooling rate per unit volume through modified 
Urea reactions is 



new 1020 rrnS I -3 -1 

— ^ ~ 9.6 X 10 Tg p(r, erg cm s 



(39) 



Equating the two rates gives the equilibrium temperature 

.,..e.ex,o-(j^)*(g)-*(^)-' K,«, 

Note that T^q is higher than Tt, if B > 7 x 
lO^'^ (pi5/0.7)"/^''(L/2km)-i/** G. Fields larger than that 



would thus be able to dissipate enough energy and balance 
neutrino cooling even in the early phase when the solenoidal 
and irrotational mode are still degenerate. This conclusion 
describes a regime that was not considered in Thompson 
& Duncan (1996): very strong magnetic fields (~ lO^'' G) 
decaying and heating a NS core at very high temperatures 
(~ 10® K) and at very young ages (years to centuries). The 
resulting evolution of magnetic field and temperature are 
coupled, as already shown by T D96. O ur solution (eq. I40|l 
joins smoothly the one found bv lTD96l (their eq. 31) in the 
sense that both give the same value of the magnetic field 
strength (B « 7 x 10^^ G) when calculated at Tt,. The two 
regimes are, in this sense, complementary, forming a con- 
tinuous evolutionary sequence through Ttr for an arbitrarily 
large magnetic field whose decay is driven by ambipolar dif- 
fusion. 

In order to better illustrate this, we calculate here the joint 
evolution of the magnetic field strength and the equilibrium 
temperature, according to the equilibrium conditions dis- 
cussed above. Consider the rate of magnetic energy dissipa- 
tion per unit volume 



B_dB 



Ant 
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(41) 



Inserting eq. (|36p and (|40|l in (|41|) gives 



dBu 
dt 



-3.12 X 10" 



whose solution is 



1.12 X 10" 



(ft)"' 



— + ( — !— 



(43) 



with Bi the strength of the magnetic field at the initial time 
ti. As an illustrative example, in Fig. |4] we show the evo- 
lution of the core magnetic field according to eq. (|43|l . for 
three different initial values. The corresponding equilibrium 
temperature evolution is: 



6.6 



1.1 X 10" 



t 



+ 



1 



B,. 



(44) 



We caution that the equilibrium temperature as a function 
of density is in principle different from the actual tempera- 
ture profile throughout the NS core. Finding this would re- 
quire solving the heat flux problem self-consistently, which 
includes account for the strong suppression of heat transport 
across magnetic fleld lines in a superstrong magnetic fleld. 

With the above expressions for the equilibrium regime, 
the decay timescales of the two modes of ambipolar diffusion 
can be evaluated self-consistently - given an initial magnetic 
core field Bt,i - and their values compared. In Fig. [5] we 
show, for illustration, the evolution of the two timescales for 
an initial (uniform) core magnetic field Bt,i = 5 x 10^^ G 
and at a given value of the density prs = 0.7. Clearly, field 
evolution is always determined by the longest decay time: 
as long as particle friction dominates, both modes decay on 
the timescale i;^J'°-^'y\ Once chemical equilibrium imbalance 
overtakes particle friction (at t « 10* yrs in our example), 
the two modes split: the solenoidal mode is unaffected and 
continues to evolve on the (now shorter) time t^J'"'^'^^ ^ while 
the irrotational mode is now subject to a slower evolution, 
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Figure 5. Evolution of the ambipolar diffusion decay timescales 
determined by each of the two friction mechanisms, separately, 
calculated at a particular density p = 2.5p„uc and for a given 
initial field strength Bi = 5 X 10^® G. The timescales evolve ac- 
cording to their dependence on the temperature T and magnetic 
field B. The evolution of T(t) and B(t) is given by eqs. I I44I I and 
1143 l l and is thus determined by np-friction alone, initially. How- 
ever, the /3-reaction timescale is a much stronger function of time 
than the np-friction timescale, in this regime. Curves are dotted 
where the two timescales become comparable, within a factor of 
a few, and the approximation under which timescales are calcu- 
lated - through eqs. 1143 l l and 1 144 l l - does not hold anymore. At 
this point, the temperature is close to the transition temperature 
Ttr and the magnetic field close to the value 7 X 10^^ G (see 
text) . 



determined by the efficiency of /3-reactions. 
We show in Fig. |4]the evolution of the core magnetic field in 
the regime described here, for three different initial values. 
More details are given in the captions. 

We can also estimate the time At after which the 
temperature, once equilibrium between heating and cooling 
holds, reaches the transition valu43 Ttr 



1.36 X 10** 



Pl5 

0.7 



f— y 



T \ — 

L \ 35 



(45) 



This gives At ~ 4500 yrs for Bi = 10^^ G and has an asymp- 
totic value (Af)„,ax — 1.36 x 10* yrs for very large B. For 
Bi =3x10^*^ G, At~ 1.1x10* yrs, ~ 80% of the asymptotic 
value. The comparison with a NS standard cooling scenario, 
without any heat sources, is striking, since in this case Tt,. 
is reached in somewhat less than 10 yrs. 
Once the field has decayed to the limiting value ~ 7 x 10^^ 
G, the temperature equals Tt, and the results obtained by 
ITD96| hold. We conclude that the early thermal evolution of 
the core can be significantly altered (slowed down) by ultra- 
strong field decay. A detailed analysis of the consequences 



This time is zero by definition for Bi = 7 X lO^^ G 
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Figure 6. Value of the equilibrium temperature (eq. I44|l as a 
function of density throughout the NS core, for four selected val- 
ues of the magnetic field strength, as indicated in the figure. 



of this conclusion is beyond our scope here and will be the 
subject of future study. 



5.1 Comparison with recent studies of core 
heating in magnetars 

Surface temperatures inferred from the X-ray spectra of 
SGRs/AXPs (> 0.35 keV, or > 4 x 10^ K) are a factor 
~ 2 3 higher than expected from cooling calculat ions of 
NSs with their estimated ages, ~ 10* yrs (Kamin ker et al.l 
I2OO7I) . The heat source internal to the NS needs thus not 
only be strong enough to provide the excess energy but it 
also needs to be effic i ently converted to surface heating. 

iKaminker et al.l (|2007 *) considered how the surface tem- 
perature of a cooling NS would be affected when a (slowly 
decaying) heating source in its interior were included, in 
order to account for the heating provided by the decaying 
magnetic field. These authors suggest that, if heat is released 
in the core or deep crustal layers, only a minor fraction of it 
is transported to the surface and radiated thereby as X- 
ray photons. Rahter, most of the released heat produces 
an enhancement in neutrino emission - which is strongly 
temperature dependent - and is thus radiated locally (via 
neutrinos). In fact, in their calculations the local temper- 
ature in the core is affected only slightly by internal heat- 
ing. Core temperatures never reach 10^ K and the surface 
can't ever be kept as warm a s AXP s/SGRs surface emis- 
sion indicates. Kaminke r et al.l l|2007h conclude that heat is 
most likely released in the outer regions of the crust, where 
neutrino-emitting processes are much less efficient and heat 
is thus most efficiently transported to the surface, with only 
minor losses. 

Their results are therefore at variance with ours, that 
envisage an efficient heating of the core up to very high 
temperatures (as a function of the magnetic field strength). 
The difference bewteen the two conclusions is readily found 
in the parametrization of the heating source chosen by 
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iKaminker et all |20o3), which corresponds to a different 
physical scenario. 

Their model heat source is independent on the mag- 
netic field strength and on the NS internal temperature. The 
maximum initial heating rate they consider (Qo = 3 x 10^'^ 
erg/s) effectively corresponds to our expression (|38p with 
B~ 1.6 X 10^^ G. The heating rate in their model decreases 
exponentially in time, with a fixed time constant compara- 
ble to, but somewhat longer than, the source estimated age. 
As opposed to this, we have i) allowed for initially stronger 
fields and, thus, larger heating rates and ii) considered a 
time-varying decay timescale, since the ambipolar diffusion- 
driven field decay is a function of the NS internal tempera- 
ture and field strength (eg. I38p . In particular, the values of 
T and B at each epoch are determined self-consistently by 
the equilibrium condition between cooling and he ating. 

As a consequence, in the model by Kaminker et al.l 

l|2007ll the released heat does not affect the value of Q while 
it enhances neutrino emission. The latter being a strong 
function of the temperature, the net result is a slight en- 
hancement of the interior temperature, to make neutrinos 
able to carry away almost all the eccess heat. 

In ambipolar-diffusion driven field decay, on the other 
hand, heating has a feedback on both cooling and heat- 
ing itself. As long as the temperature is very high the NS 
cools, with heating providing just a minor perturbation to 
the dominant process of z^-cooling. As the core temperature 
drops, however, the field decay rate grows (eq. I38p while 
neutrino emission drops. Eventually, the two rates become 
almost equal and at th is stag e their equilibration plays a 
key role, as stressed bv ITDQSI . Near equilibrium, tempera- 
ture variations have a strong feedback on both heating and 
cooling - and with opposite effects. This forces the tempera- 
ture toward T^, (eg. l4Up : as the field dissipates the temper- 
ature drops slighlty and a new equilibrium between heating 
and cooling is reached, at a slightly smaller temperature and 
field strength. The equilibrium temperature within the NS 
core, as a function of density, obtained through our eq. (|40p 
is shown in Fig. (|6]) for four different values of the average 
core magnetic field strength. In Fig. ((T)) we show, for a given 
initial magnetic field strength, the equilibrium temperature 
throughout the magnetar core at four different epochs of its 
evolution according to eq. (|4ip . More details are given in the 
caption. 

Therefore, the internal heating source considered by 
IKaminker et al.l (|2007l ) differs from ambipolar diffusion- 
driven field decay. When ambipolar diffusion in the core - 
and the associated heating - is taken into account, the NS 
core can remain at fairly high temperatures for a long time 
if the decaying field is around ~ 10^® G. 



6 CONCLUSIONS 

In this paper we have investigated some implications of one 
of the key ansatz of the magnetar model; namely, that mag- 
netars do form with millisecond spin periods and a (mainly) 
toroidal magnetic field, generated through the strong differ- 
enti al rota tion of the collapsing proto-rieutron star (|DT92| 
and lTD93l ). 

Building on our earlier work l|Stella et al. I l2005l . 
iDaU'Osso fc Stellall200'iD , we showed that one major imph- 
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Figure 7. The evolution of the temperature profile within a mag- 
netar core (10^'' g cm~'^ < P ^ 10^^ g cm~^) for a specific initial 
value of the toroidal field strength (Bt = 4 X 10^^ G), corre- 
sponding to a total magnetic energy Eb — 4.6 X 10^" erg or an 
ellipticity cb ~ —4.6 X 10~^. The four curves describe the temper- 
ature profile at four different epochs, the thick solid curve at the 
bottom defines the temperature at (and below) which chemical 
equilibrium imbalance becomes the major limiting factor for am- 
bipolar diffusion of the irrotational mode (the regime described in 
TD96). Above the thick line, the irrotational and solenoidal mode 
are degenerate and the regime described in the previous section 
holds. 



cation of this scenario is that such objects can become strong 
sources of GWs in the first few days after formation. This 
results if there exists at least a tiny misalignment (angle 
x) between the rotation axis and the symmetry axis of the 
magnetic field, at birth. The newly born NS is distorted to 
a prolate shape by the toroidal field, and is freely precess- 
ing because of x 7^ 0. Under these circumstances internal 
viscous dissipation of the precessional motion will drive the 
magnetic symmetry axis orthogonal to the spin axis, the 
most favourable geometry for GW emission. 

We discussed current uncertainties in various aspects of 
the model and introduced simple approximations to treat 
each of them. We developed a simple analytical model de- 
scribing the early rotational evolution of newly formed mag- 
netars, that includes an ideal magnetic dipole torque plus 
the GW torque acting on an orthogonal, prolate rotator (see 
eq. |4]). We then estimated the magnitude of GW emission 
from magnetars s as a function of their initial spin, internal 
magnetic energy and external (dipole) magnetic field (as well 
as their mass and radius). 

Our main conclusions can summarized as follows: 

• if magnetars are born with spin period less than 3 ms, 
internal toroidal fields ^ 3 x 10^^ G and external dipole fields 

2 X 10^* G then the expected GW signal would be strong 
enough to be detectable with Advanced LIGO/ Virgo class 
detectors out to the Virgo cluster, where their formation 
rate may be ~ 1 per year; 

• the estimated optimal S /N ratios for match-filtered sig- 
nal searches (with one detector only) are very encouraging. 
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However, as already noted by IStella et alj (|2005l ). optimal 
signal searches have unafFordable computational costs. The 
developement of sub-optimal signal search strategies is re- 
quired, a task that is currently under way. 

• if our scenario holds, the rotational energy ~ 3 x 10^^ 
(P/ms)^ ergs of the newly formed magnetar will be emit- 
ted mostly as GWs, in the first few days after formation. 
As a consequence. Supernova Remnants (SNRs) surround- 
ing evolved magnetars would not be expected to bear the 
signature of an excess energy injection (> 10'''^ erg) soon 
after formation, since GWs do not interact appreciably with 
the expanding shell. 

In particular, the condition that a newly formed mag- 
netar be detectable as a GW source from Virgo-cluster dis- 
tances turns out to be almost equivalent to the condition 
that it radiate less than lO^'^ erg through magnetic dipole 
radiation. The two requirements are met for nearly over- 
lapping regions of the internal and external magnetic fields 
parameter space. Stat ed differently, if t he Ga lactic magnetar 
candidates studied bv lVink fc Kuipej ()2006l l were born with 
millisecond spin periods, Bt > 3 x 10^*^ G and Bd ^ 2 x lO", 
they would have lost most of their rotational energy through 
GW emission. Their SNRs would thus not show any excess 
energy compared to other SNRs - as observed, indeed - and 
the GW signal they emitted could have been detected out to 
~ 20 Mpc with Advanced LIGO/Virgo class interferometers. 

Finally, we considered the evolution of an internal field 
> lO^'' G as a r esult of am bipolar diffusion, as already envis- 
aged bv lGR92l and lTPgd Our aim here was twofold: first, 
we investiga te this high B-field regime, for which the cal- 
culations by ITD96| are not appropriate. Second, we showed 
that even fields this strong have (at least) a slow decay mode 
through ambipolar diffusion, that is active soon after forma- 
tion. This process can prevent the cooling of the magnetar 
core below a temperature of ~ 10® K for hundreds to thou- 
sands years. This conclusion is expected to have significant 
implications for our understanding of AXPs/SGRs. 



APPENDIX A: SPIN ENERGY AND 
FREEBODY PRECESSION ENERGY OF A 
ROTATING ELLIPSOID 

We give here a quick derivation of the expression for the 
energy of freebody precession of a fluid star subject to 
both centrifugal and magnetic deformations. A more gen- 
er al discussion , in the conte xt of "ordinary" NSs, is found 
in I Jones fc Anderssonl l|200lh . 

First of all, we write down the frequency of the fr eebody 
precession mode derived bv lMestel fc Takharl l|l972f ) 

h - h 

Wpro = ——: flcosip = esficost/i , (Al) 

h 

where angles throughout this section are those defined in 

Fig.EB 

Following I Jones fc Anderssonl[200ll (and references therein), 
we define the moment of inertia tensor of the fiuid NS as a 
linear combination of three contributions, namely a (spheri- 
cal) gravitational part, plus two axisimmetric perturbations 
provided, in our case, by the centrifugal and magnetic fields, 
respectively. Hence: 



L 




Figure Al. Schematic representation of the angles of interest to 
the problem. L represents the invariant angular momentum vec- 
tor, Q is the angular frequency vector and n^ is represents the 
unit vector along the simm etry axis of the magnetic fi eld. The an- 
gle (/> ~ egsin-)/)cosi/i <^ i/i l lJones fc Anderssm]|200ll) . Therefore, 

I = IoS + AIninnnn - S/3) + AlBingnB ~ 5/3) (A2) 

where d_ is the identity, nn,s represent the unit vectors along 
the spin and magnetic field axis and AIq^b are the magni- 
tudes of the corresponding perturbatio ns of the i nertia ten- 
sor. The eigenvalues of the latter are l|Jones fc Anderssonl 
iool): 

Ii = lo + -Ain - -Alfl 

2 a 1a 
h = lo + -Ain - -Alfl 

Is = Io + |AlQ-f |Alfl =Ii-|-AIb (A3) 

These expressions show that the centrifugal deformation 
modifies all eigenvalues in the same way, as if it was effec- 
tively a spherically symmetric, additive term. The magnetic 
deformation, on the other hand, does introduce an asim- 
metry between the first two eigenvalues and the third one, 
inducing the magnetic ellipticity es = (I3 — Ii)/li. The an- 
gle (j) between CI and the angular momentum axis is always 
smaller than the tilt angle x of ths magnetic symmetry axis 
(iJones fc Anderssonir200lD . In formulae, to first order in es 
it can be shown that (by using eq. IA2|I 

sint/) ~ ~ efl sin?/'cos'!/' ^ , (A4) 

from which 'ip — x ~ — X- 

Following the argument bv ICutler fc Jones! (120011 ). the NS 
angular momentum can be written as i = Ifl and its kinetic 
energy as = (l/2)/f2 ■ ft, from which (to first order in 

L = (1 -|-2eBCosV)^''^ ^ Ii^^ (1 + 2escos\)^''^ 

Ek = ili!^2(l-FescosV) ~ ilin^(l + escos\). (A5) 
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From the latter equation we can obtain the energy of free- 
body precession by taking the difference between the sec- 
ond equation in (|A5|l and the spin energy of the same ellip- 
soid that, with the same total angular momentum L, spins 
around the axis having the greatest moment of inertia. The 
latter configuration is indeed the one that minimizes the en- 
ergy, at constant angular momentum. It is thus the one to- 
wards which a freely precessi ng spheroid will evolve, g i ven an 
internal dissipatiy e process (jMestel fc Takharl Il972l : Ijoned 
ll976l : ICutleill2002l) . Therefore, if Q is the angular frequency 
of the precessing spheroid and Qf the angular frequency once 
freebody precession is completely damped, conservation of 
angular momentum implies 



(A6) 



Writing the minimum spin energy as E„i„ = ^Imax^Fi the 
freebody precession energy is 

1 2 

Eprcc = Ek — Imax^^F 

= ^l^n'' ^^l + eBCOs\- :^{l + 2eB)Y (A7) 

This general expression can be specialized to the case of an 
oblate (Iniax = I3) or a prolate (l„,ax = Ii = I2) ellipsoid 
giving^, respectively 

Eprcc — — Ilf^^es sin^;^ oblate ellipsoid (eg > 0) 
Ep,cc — ~2^l^^^s cos^;^ prolate ellipsoid (eg < 0). (A8) 
to first order in e_B . 

Self consistence of the above is warranted by Epr^c — * for 
X — > in the oblate case, and for x ^ ''r/2 in the prolate 
case. 

From the above formulae we can eventually relate the time 
derivative of the freebody precession energy to the time 
derivative of the tilt angle x, in the case of a prolate ellip- 
soid. Since conservation of angular momentum is required, 
Q changes only as a consequence of changes in x (cfr. eq. 
IA6[I . which makes the time derivative of Ep^cc a function of 
X only. Taking the time derivatives of the first of eq. HA5|) 
and the second of eq. (|A8|I and requiring angular momentum 
conservation we obtain (to first order in es) 

dEp.ec 



dt 



(A9) 



where was defined in § 12.1.41 



APPENDIX B: CALCULATION OF THE 
DAMPING TIME OF FREEBODY PRECESSION 
THROUGH BULK VISCOSITY 

In this appendix we describe the calculation that leads to the 
estimated timescale for dissipation of the free precessional 
motion through bulk viscosity (eq. I13p . According to the 
definition of Td (eq. I12|l . we need an expression for both 
the bulk viscosity coefficient, (, and the precession-induced 
density perturbation, Sp. 

Note that the preeession energy is positive in both cases, as it 
must be 



For the bulk viscosity coefficient, we recall here the gen- 
eral expression (fTTjl 

nr(dp/dx)^dx/dn n{dp/dx)^dx/dn 

t<.e(C,J = ~ 5 . {HI} 

1 -I- (tjrj^ LO'^T 

For npe matter, as we have assumed through- 
out, the timescale for /3-reactions (t/j) is (cfr. 
iReisenegger fc Goldreichlll992l ) 



rp = 



3n„ 



0.23 



(B2) 



Here x is the equilibrium fraction of charged particles (see 
below), Ef^ = (/iV2?n„)(37r^n„)^/^ is the neutron Fermi 
energy and ~ 5 x lO^^Tg (p//9„ac)^''^ erg"^ s~^ is the rate 
of /3-reactions when the combined fluid is out of chemical 
equilibrium by an amount 5p, = Spn — 5 fie — Sfip. Finally, 
p„uc ~ 2.8 X 10^* g cm"'' is the nuclear saturation density. 
The last step in eq. holds for ujpiaT/s S> 1, the relevant 
approximation here. 

The pressure of degenerate npe matter, neglecting the small 
contribution from protons, is given by p ~ |m„_EF„ + 
jTrie-EFc where the equilibrium proton (and electron) frac- 
tion is (ReisoncEEcr & Goldrcich 1992) 



/np_\ 



!^^6xl0-^^ 



(B3) 



From all the above, one obtains the required expression (cfr. 
ISawycr 989): 



6 x 10 



-59 P J 



cos^x VO.' 



As a starting poin to determine the density perturbat i on Sp , 
we refer to the expression given bv lMestel fc Takhaij (Il972l ) 
of the centrifugal distortion of a fluid star, pn, as a function 
of a spherical coordinate system whose origin is at the star 
center and whose pole is the magnetic pole (r, 6, A) . 
The amplitude of internal field of motion, 4, is determined 
by the non-spherical part of the rotational distortion of the 
fiuid. The perturbation of the density profile is 

5pn = /5n(A-nt)-pn(A) = i/(r)A'(x,e,A,nt) . (B5) 

where the function K(x, 9, A, Q.) is 

k(x,e,\,n) = sin^x[i - ^'2(m)] X 

X {sin2A sin2nt - cos2A(1 - cos2f2t)} -I- 

-I- 3 sinxcosxsin2S [sinA sinf2t — 

- cosA(l - cosM)] . (B6) 

Here, P2(p) ~ P2{cos6) is the Legendre polynomial. The 
radial function /(r), that depends on the stellar model as- 
sumed, can only be obtained analytically through approxi- 
mate (and/or idealized) calculations. 

In particular, the centrifugal deformation of a fluid mass 
with a polytropic equation of state is an old problem that 
has been largel y discussed in the liter ature. We refer, here, 
to the study by IChandrasekhail (|l933l ). where both analyti- 
cal formulae and numerical values for the relevant functions 
are derived. 
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Given the polytropic EOS p = 
adimensional density variable 



en 



^pi+i/n^ one can define an 



(B7) 



The ratio of the rotation energy to the gravitational binding 
energy defines the rotation parameter v 

v = ^^ . (B8) 
2-KGpc 

in powers of which the adimensional density O can be ex- 
panded. In the above G is the gravitational constant and pc 
the central density of the fluid star. For w <C 1, a first order 
expansion will sufBce. Introducing the adimensional radial 
coordinate ^, such that r — a^, the first order expansion of 
the density profile becomes 



e(c) = eo(e)+«^(e) 



(B9) 



where 80 is the unperturbed density profile, solution to the 
Lane-Emden equation for a non-rotating polytrope. The pa- 
rameter a above is related to the polytropic equation of state 

^ \2tvGJ 



jn + l) k 



4tvG 



(BIO) 



The problem is thus reduced to finding an appropriate ex- 
pression for the function V'(C)- -fri particular, for polytropes 
of index n the following expansion holds 



9(0 ^eo + v [MO + A2MOP2M] 



(Bll) 



where ipoiO the spherical part of the centrifugal defor- 
mation and the non-spherical part of the deformation cor- 
responds to the second term in square parenthesis. The co- 
efficient A2 is 



A2 



en 



, , (B12) 

where a£_R = R and the prime denotes a first derivative. 
An approximate analytical expression for t p2(0 is given as 
an ei ghth-order polynomial expansion in ^ (|Chandrasekhail 
1 19331 ) 



14^ 504^ 



(B13) 



33264 

which formally completes the problem at hand. 
As it has been shown by numerical integration of the equilib- 
rium equations, a first-order expansion in v provides a rela- 
tively good approximation to slowly rotating polytropes; in 
the case n = 1, v ^ 0.075 gives approximately the rotation 
rate at which deviations from the first ord er perturbation 
theory become non-negligible (|Tassoullll978l ). In particular, 
direct integration of the equilibrium equations reveals a sys- 
tematically larger deformation of rapid rotators compared 
to the results obtained through the first order perturbative 
approximation. Therefore, our general expectation would 
be that the latter underestimates the real deformation of 
a rapidly rotating NS {5p), thus underestimating the effi- 
ciency of bulk viscous dissipation. 

The polytropic index n = 1 greatly simplifies the analytical 
treatment and allows one to derive quite straightforward for- 
mulae. The adimensional radius of the polytropic star will 
be = TT and the nur nerical values oi i l}2(£.R ) and i/'2(Cfl) 
have been tabulated bv lChandrasekhaJ (119331 ). allowing to 



derive A2{n = 1) ~ —0.54833. Finally, a is determined by 
the linear scale for the NS radius so that: 



a ~ 3.82 X 10' 



V12 Km/ 



(B14) 



and the EOS thus becomes p = 6.1146 x 10* (R/12 Km)^ . 
The density perturbation of eq. (|B5|l . the one that deter- 
mines the precession-induced internal motions damped by 
bulk viscosity, amounts to only the non-spherical part of the 
centrifugal deformation. The spherical part does not cause 
any radial periodic motions and, thus, does not perturb the 
local chemical potential equilibrium of NS matter. We are 
therefore left with: 



(B15) 



where the angular part corresponds to the function 
K(x,0,\,Q,). This must be inserted in eq. (|B5|l to obtain 
the full expression for the non-spherical density perturba- 
tion. Recalling eq. (|12|l . we thus get the full expression for 
the bulk viscosity damping rate: 



EdiBs = OOTio Lo^ 



where: 



dV 



5p 



p I LJ^ 



= eOTi^o / dV(<5p)^(B16) 



{p.vfAl^l{i)K\x,e:\^) 



(B17) 



so that, eventually, the full expression for Edisa can be writ- 
ten as: 



Edisa 



2.6 X 10"rfon**a'' / V2(C)rdC / ^ sinSdSdA 

3 



2.2xl0« f^y/msy/^y 



X / '/'2(C)rdC / A' sin^dedA 



(B18) 



The radial integral, solved from ^ = to = tt gives 138.04, 
while the angular part gives 

(247r/5)sin2x(l + 3cos2x)- Even- 
tually: 

Ed.aa ~ 4.7 X lO"*** sin'x(l + 3cos'x) 

(t5?^)"(t)'(uL)'--- '™) 

As a last step, given the energy of the freebody precession 
mode derived in § \^ we obtain the dissipation timescale 
(see eq. [13}: 

13.5 cos^X / Eb a V P a ^ 
UOSOerg/ \ms) 



2£p 



Td 



Edisa 



sin^x(l + icos^x) 
( T ( M 



V 1010k 



\AMr. 
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